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~''';^Thi8  report  covers  some  aspects  of  the  formulation  of  the  stabilization 
problem  by  aerodynamic  means.  Results  of  experimental  work  are  not  included  but 
are  referred  to  in  appropriate  places.  The  classical  solution  for  ambipolar 
diffusion  in  tne  positive  column  of  the  glow  discharge  must  be  extended  to  re- 
flect the  influence  of  current  and  pressure.  The  influence  of  current  is  given 
by  the  inclusion  of  appropriate  energy  equations  and  the  pressure  influence  by 
tne  addition  of  two-body  recombination.  It  was  found  that,  from  dimensional 
analysis,  the  characteristic  lengths  for  ambipolar  and  heat  diffusion  for  
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I.  INTRODUCTION 


In  the  high  energy  laser  field,  the  electric  discharge  convection  laser 
(EDCL)  offers  good  efficiency  as  well  as  simplicity  of  handling  relative  to 
gasdynamic  and  chemical  lasers.  Progress  towards  higher  laser  outputs  from 
EOCLa,  however,  has  been  Impeded  by  the  collapse  of  the  glow  discharge  at  the 
Increased  currents  and  pressures  desired.  This  glow  collapse  (usually  re- 
ferred to  as  arcing  or  constriction)  represents  the  practical  upper  limit  of 

the  electrical  power  that  may  be  pumped  Into  the  necessary  non-equlllbrlum 
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processes.  It  has  been  shown  that  aerodynamic  stabilization  can  be  used 
very  successfully  In  overcoming  the  glow-collapse  limits.  It  should  be 
pointed  out  here  that  the  laser  application  Is  but  one  application  of  the 
relatively  new  Interdisciplinary  field  of  cold  plasma  chemistry. 

Aerodynamic  stabilization  Is  based  on  two  separate  effects  usually  re- 
ferred to  as  convection  (laminar)  and  mixing.  The  purely  convective  effect 
Is  said  to  decrease  the  residence  time  of  the  particles  through  the  discharge 
and  thereby  to  counteract  unstable  mechanism  that  lead  to  the  glow  collapse^. 
The  mixing  effect,  which  may  be  due  to  vortex  flows  or  to  turbulence  from 
boundary  layers  or  grids,  has  been  observed  to  have  powerful  stabilizing 
characteristics.  Convection  scales  up  with  flow  velocity  whereas  mixing 
depends  on  the  Intensity,  spectral  distribution,  and  location  of  the  eddies 
or  vortices. 
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II.  SPECIFICS 


Ac  Che  Technlon,  che  auchor  had  Che  opporcunlcy  co  parclclpace  in 

experlmencal  work  dealing  wlch  vorcex-flow  sCablllzaclon^*^  and  wlch  super- 
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sonic  flow  sCablllzaClon  ’ . These  Cechnlques  complemenC  che  work  ChaC  has 
been  underway  aC  che  Poscgraduace  School  on  turbulence  stabilization  of 
discharges.  Some  results  of  Che  Technlon  work  have  already  appeared  (see  Refs. 
6-11)  and  more  will  be  forthcoming. 

A significant  portion  of  Che  author's  time  was  dedicated  to  a study  of 
theoretical  aspects  of  discharge  stabilization.  The  author  was  privileged  to 
colaborate  with  Dr.  Y.  Khalt  In  establishing  a perspective  regarding  some  of 
the  author's  previous  results  with  turbulence  stabilization.  This  effort 
resulted  In  Che  writing  of  a paper  titled  "Influence  of  Turbulence  on  Diffuse 
Electrical  Gas  Discharges  Under  Moderate  Pressures,"  which  Is  being  submitted 
for  publication.  Ocher  theoretical  work  Involved  a close  look  at  Che 
SchoCcky  solution  for  che  positive  column  In  a glow  discharge  with  the  aim 
at  establishing  Its  use  for  a stability  formulation.  This  effort  was  fruit- 
ful on  two  Intermediate  aspects  of  the  problem  and  these  are  Included  here 
as  Appendixes  A and  B. 

A proposal  was  written  to  che  United  States-Israel  Blnatlonal  Science 
Foundation^^  which  extensively  outlines  the  progress  of  the  entire  effort. 

If  awarded,  the  cooperation  of  the  Technlon  and  Che  Postgraduate  School 
would  be  formalized  (however,  the  budget  would  provide  neither  salaries  for 
principal  Investigators  nor  equipment  for  the  USA  experiments) . 
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III.  DISCUSSION 

! 

The  classical  solution  to  the  amblpolar  diffusion  of  charges  In  the  , 

12  13  - 

positive  column  * of  a glow  discharge  offers  no  hints  as  to  how  strlatlons  | 

and  constriction,  among  other  Important  observations,  manifest  themselves  In  I 

the  discharge.  As  pressure  and/or  current  ate  increased,  Instabilities  grow 

In  the  discharge  and  It  Is  Important  to  have  a sufficiently  complete  model 

of  the  physics  which  reflects  these  effects.  In  other  words,  before  flow 

effects  such  as  (laminar)  convection  and  turbulence  are  Introduced  It  Is 

necessary  to  re-examine  some  assumptions  Inherent  In  the  Schottky  solution. 

A partial  list  of  phenomena  not  reflected  In  the  Schottky  solution  may 

be  given  as  follows: 

1}  effect  of  sheaths  at  the  walls 

2)  effect  of  electrodes  and  Intermediate  regions 

3)  effect  of  recombination  In  the  gas  (pressure  and  composition) 

4)  effects  of  the  energy  balance.  Including  radiation  and  the  elevation 
of  the  electron  temperature 

The  list  can  be  rather  long.  We  closely  examined  the  energy  balance  because 
the  Instabilities  are  believed  to  be  thermal  disturbances  (l.e.,  the  tempera- 
ture becomes  non-homogeneous)  and  recombination  In  the  gas  because  It  closely 
reflects  the  effect  of  pressure  and  composition.  For  simplicity,  a two- 
dimensional  (Infinite  parallel  plate)  geometry  was  considered,  see  Fig.  1. 

When  is  constant,  7*J  ■ 0 Implies  that  ■ J(y)  or  that 

n ■ n(y)  because  ■ 0 ; moreover,  V x E ■ 0 Implies  that  E^  ■ E(y)  . 

The  quaslneutral  equations  for  steady  state  may  be  written  as^^*^^ 
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Figure  1.  Geometry  of  Positive  Column 


■ [2II7  “ ''l^  “ 

(1) 

- 

- 3/3y  [K  3T/3y]  ■ J*E  ■ enu^E^ 

(2) 

• 

2 — V 1 k(T  - T)  - eu_E,^ 

m^  evi  4 B 

p - NkT 

(4) 

with  boundary  conditions 

at  y ■ 0 3n/3y  ■ 3T/3y  ■ 0 

at  y*+a  n*0  and  T ■ T 

-•  o 

In  the  above  formulation  x Is  the  axial  coordinate  and  y the  Interplate 
coordinate,  n the  charge  concentration  In  the  amblpolar  region,  T the 
electron  temperature  and  T the  gas/lon  temperature,  K Is  the  concentra- 
tion of  the  neutrals,  D the  coefficient  of  amblpolar  diffusion,  K the 
thermal  conductivity,  m^  the  electron  mobility,  and  a the  Ionization 

and  recombination  coefficients  respectively,  E the  axial  electric  field, 
the  collision  frequency  between  electrons  and  neutrals,  e the  elec- 
tronic charge,  p the  pressure,  and  m^  and  m^  the  mass  of  the  electrons 
and  neutral  respectively.  There  are  several  assumptions  Inherent  In  this 
formulation,  namely  that  N>>n  , T^>T  , and  that  p ■ constant  . In  addition, 

we  are  using  a symmetry  condition  In  the  boundary  conditions  and  we  are 

12  13 

neglecting  of  wall  space  charges.  Glow  discharge  data  * , In  addition. 
Indicate  that  as  well  as  appear  to  be  constant  In  the  discharge  at 

least  before  constricting.  The  electron  energy  equation,  Eqn.  3,  should  have  to 
reflect  the  energy  dependence  of  the  collision  cross  sections  and.  In  addition. 
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noa-elasclc  collisions  — since  we  are  Interested  In  laser  operation. 
However,  we  shall  not  deal  with  Eqns.  (3)  and  (A). 

Using  simple  kinetic  theory  to  describe  the  gas  temperature  'ependance 
of  the  coefficients  (D  , u , K , etc.)  and  defining  a stationary  state 
where 


or 


v^n*  “ an* 


n*  “ v^/a 


(5) 


we  can  proceed  to  nondlmenslonallze  Eqns.  (1)  and  (2).  We  shall,  however, 
not  use  the  plate  spacing  to  render  y dimensionless,  but  Introduce  two 

The  gas  temperature  is  made  non- 


characteristic lengths,  and 

dimensional  by  the  use  of  a wall  temperature  T 


are: 


- d/dy  ^)  - n - n^ 


The  resulting  equations 


(6) 


Ln  = 


T D 
e ao 

2v.  T 
lo  o 


d/dy  (T^''^  1^)  - nT 


^T  ~ 


T K 
o o 


(7) 


(8) 


(9) 


It  seems  reasonable  to  expect  that  In  a conventional  (unconstricted) 
glow,  the  two  characteristic  lengths  be  about  equal  to  each  other  and  that 

L /a  K 0(1)  and  L_/a  a 0(1)  (10) 

n 1 


This  expectation  is  also  required  by  the  eigenvalue  nature  of  Eqns.  (6)  and 
(8).  It  may  be  shovm  (see  Appendix  6),  that  except  for  the  equilibrium  limit, 
Eqn.  (6)  Is  an  eigenvalue  equation,  (the  classical  solution  being  a cosine 
when  no  recombination  is  present  and  when  T % constant) . me  nature  of 
Eqn.  (6)  Is  discussed  further  at  the  end  of  this  section. 

In  the  light  of  Eqn.  (10) , It  Is  Interesting  to  look  at  some  actual 
values  of  and  . Particularly  extensive  Information  for  helium  Is 

reported  In  Ref.  16,  and  this  Is  used  In  estimating  the  values  shown  In 
Fig.  2,  which  are  limited  to  the  range  of  applicability  of  the  equations 
In  Ref.  16.  An  appropriate  form  for  the  coefficients  Is  obtained  by  Intro- 
ducing Che  gas  density  , 


L 

n 


N 


0 


V 2 <.\/V 


‘•T»o 


(E/N  ) 
o 


-1 


T K a 
o o 


eo  o 10 


(11) 


(12) 


i 


The  above  forms  show  "similarity"  since  LN  is  a function  of  E/N  only. 

Examining  Fig.  2 , It  Is  evident  that  L^N  Is  not  equal  to  L^N  , certainly 

not  for  large  E/N  , wheareas  follows  closely  data  reported  for  helium 
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In  Che  literature  * (the  fact  that  we  are  looking  at  parallel  plates  and 
Che  data  Is  for  Cubes  Is  Immaterial) . Since  In  a discharge  tube  n and  T 
share  a common  diameter.  It  Is  necessary  to  reconcile  the  discrepancy  of  the 
two  curves. 

14 

It  Is  possible,  chat  the  assumption  which  brings  T Into  Eqn.  (1) 

and,  therefore,  into  Eqn  (11)  la  Invalid.  But,  even  If  v-e  were  to  eliminate 
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the  factor  from  the  equation,  the  curves  would  diverge  at  the  higher 

values  of  E/N  . Constriction  of  the  positive  column  Is  another  possibility 

but  in  this  case  L-  would  be  expected  to  be  greater  than  L for  any  given 
T n 

E/N  and  just  the  opposite  seems  to  be  true.  Uncertainties  in  the  estimate 
of  a , the  recombination  coefficient,  would  affect  the  value  of  but 

this  estimate  is  the  one  that  follows  nicely  the  experimental  values. 

The  method  of  reconciliation  of  the  two  lengths  proposed  here  is  based 
on  the  alternative  that  the  equations  previously  written  are  incomplete  and 
that,  in  particular,  in  Eqns.  (1)  and  (7)  should  be  replaced  by  a 

parameter  (C^)  which  makes  up  for  the  difference  in  the  two  curves.  This 
matter  is  discussed  in  Appendix  A of  this  report. 

Consider  now  Eqn.  (1)  with  constant  coefficients.  We  proceed  here  to 
nondimens lonalize  it  in  a conventional  way 

- D n”  - V.  n - an^  (13) 

a 1 

n'(o)  ■ 0 and  n(j2  a)  ■ 0 


Let  n(o)  • n and  e ■ an  /v 
O 0 1 

Hence  we  get 

- n"  ■ n - en^ 

n'  (o)  ■ 0 and  n(o)  ■ 1 
n(+  a/  y ■ 0 


(14) 


Now  the  Shottky  solution  is  one  for  which  e - 0 , namely, 

- n”  ■ n 

and  n • OgCos  (j  y/a) 


(15) 

(16) 
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It  can  be  seen  therefore  that  the  characteristic  length  L should  be 


L 


(17) 


On  the  other  hand,  for  the  equilibrium  limit  e ■ 1 and  the  equation 


n 


“2 

n 


(18) 


is  of  the  "Boundary  Layer"  type,  l.e.,  it  is  not  an  eigenvalue  problem,  see 
Ref.  14  pp  149-51  for  the  equivalent  problem  with  three-body  recombination. 
The  solution  to  Eqn.  13  in  an  implicit  form  is 


y 


log  r^-iT 


ft  + 


I ;3  1,1/2 


1.317 


(19) 


or  explicitly 


“ “ 2 % 


1 _ (y  + 1*317) 


12 


1 + e 


- (y  + 1.317) 


(20) 


which  is  Indeed  of  the  boundary-layer  type.  The  intermediate  cases, 

0 < e < 1 , are  represented  by  elliptic  functions,  and  this  matter  is  dis 
cussed  in  Appendix  B of  this  report. 
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IV.  CONCLUSIONS 

Some  deficiencies  of  the  Schottky  solution  have  been  re-examined  with 
the  Intent  of  establishing  a more  useful  solution  for  the  description  of 
the  transition  to  the  glow  collapse.  At  least  two  major,  recent  attempts  In 
the  literature  (Refs.  16  and  17)  base  the  steady  solution  for  the  charge 
profile  on  the  Schottky  result  and  these  attempts  have  not  succeeded  In  pre- 
senting a complete  formulation  of  the  problem.  The  Technlon  approach^^  Is 
to  base  the  perturbation  on  a physically  complete  solution  of  the  amblpolar 
charge  profile.  This  profile  must  Include  the  Ingredients  of  constriction 
and  of  the  elevation  of  the  electron  temperature. 

As  shown  In  Appendix  A,  the  characteristic  length  scales  for  amblpolar 

diffusion  and  heat  diffusion  (conduction  only)  may  be  reconciled  If  In 

Eqn.  (7)  Is  replaced  by  where  C.,  > v.  and  C.  Is  a function  of 

^ 4 lo  z 

E/N  (l.e. , the  dependence  given  by  the  difference  of  the  two  curves  In  Fig. 
2).  What  appeared  to  be  a bonus,  namely,  a Insight  Into  the  mathematical 
description  of  strlatlons,  needs  further  study  In  order  for  It  to  make  more 
physical  sense.  The  eigenvalue  nature  of  the  diffusion  equation  Is  not 
charged  by  the  recombination  term  except  at  the  limit  when  e ■ 1.0. 
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APPENDIX  A - PROBLEM  REFORMULATION 

Equation  (1)  In  the  text  presupposes  steady  flow  along  with  n > n(y)  . 
Let  us  relax  these  suppositions  and  go  back  to  the  species  flux  equations 
which  may  be  written  aa 


II  + 7.  [n(-  E - ^ 7 n)]  - i 


II+  7.  [n(y,  E - -^7  n)]  - h 
at  in 


(Al) 

(A2) 


Combining  and  eliminating  the  electric  field  we  get  the  amblpolar  diffusion 
equation  In  the  usual  manner  which  for  constant  coefficients  is 

(A3) 


3n/3t  - D 7‘’n  - n 
a 


or  3n/3t  - D ( — y + y 
^ \3x^  3y2/ 


n - an 


(A4) 


The  above  assumes  two-body  recombination  as  well.  Now,  using  Eqn.  (A3)  to 


eliminate  7^n  from  Eqn.  (A2)  yields 


9«  , “1  '•  . 

+ ■ • n 


3t  1 - D./D 

1 a 


(A5) 


But  since  7«E  ■ 0 In  Che  amblpolar  region 


3n  ^1  r _ • 

+ n Jrs  E*7  n “ n 


3t  1 - D /D 
1 a 


(A6) 


or  3n/3t  + 


3n 


- D^/dJ  \h  3x  ^ 3y 


CiCL  V 2 

■ V,  n - an 


(A7) 
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Let  us  work  with  Eqna.  (A4)  and  (A7)  In  linear  form,  i.e.,  neglecting 
recombination.  If 

n(x.y.t)  ■ f(y)  g(x,t)  (A8) 


Then,  substituting  into  Eqn.  (A4) 

f gt  “ ^a  (fgj^  + gf")  - v^fg  - ... 


(A9) 


where 


d^f 


f"  - — I g^  - — I and  g^  - 3g/3t 

dy  3x 

Dividing  by  fg  both  sides  of  (A9) 

l/g  (8,  - ■>,  8„)  - ^ f ■ -i 

At  a given  E/N  , is  constant  so  we  may  presume  that 
l/g  (g^  - Dg  8^)  ■ constant  (or  - 
D 


(AlO) 


and 


- — f"  - constant  (or  C2) 


with  ” ^ ^2  " '^i 


(All) 


(A12) 


(A13) 


The  constants  Cj^  and  C2  are  positive  and  defined  so  that  C2  ^ 
as  required  by  Eqn.  (A12) , which  now  replaces  Eqn.  (6)  since 


- D f"  - f C„ 

a 2 


(A14) 


and 


Lf  i 
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Equation  (A14)  Is  the  new  Schottlcy  equation.  The  constant  appears  in 

the  coiq>anion  partial  differential  equation 


«t  '=1  * ■ 


Now  look  at  Eq.  (A7) . At  the  centerline  3n/3y  - 0 and  f(o)  ■ f 


(A15) 


so  that 


at  y • 0 


1 - D,/D 
^ i a, 


■ ''l  ” 


(A16) 


Introducing  n ■ fg  we  obtain 


*t  ^o  yi  - Dj^/D^y  \ *x  ^o  “ ^i  ^0* 


(A17) 


We  can  cancel  f^  in  the  linear  equation  above  and  obtain  an  equation 
in  which  no  function  of  y appears 


®t  yi  - D^/D^/  ^ " ''i® 


(A18) 


At  this  point  it  is  worthvrtiile  to  eliminate  g^  between  Eqns.  (A18)  and 


(A13)  and  obtain 


°a  *xx  ■ ^1  ® “ '^i*  ' yl  - 


* (1  - D./D  ) ■ ('■l  * ‘=1>* 

1 a 


«xx  [d^  - dJ  8x  ‘ « 


(A19) 


Equation  (A19)  can  be  solved  with 


g 


T(t)  e 


ax 


(A20) 


where  a Is  found  from 


a 


+ 4C2/D 


a 


2 


(A21) 


Note  that  D.  < D and  that  C Is  anyway  a positive  parameter. 

X A 

Equation  (A21)  Indicates  that  the  roots  are  real  and  there  Is  no  periodic 
behavior  In  x — that  strlatlons  are  not  to  be  found  with  this  description. 

In  order  for  n to  oscillate  In  x , we  must  have  a situation  In  which 
new  constants  and  C2*  are  Introduced  such  that 

Cl*  - C2*  - (A22) 

*t  -■>.*»-  <=i*  * 

and  - D f"  - - C-*  f (A24) 

a L 


here,  Eqn.  (A24)  will  result  In  (non-per Iodic)  hyperbolic  functions  In  y and 
Eqn.  (A24)  In  oscillatory  functions  In  x provided  that 


4C2* 


(A25) 


As  before  we  would  like  to  have 


h*  ^ ''i  ^^2*  " 

so  as  to  match  L to  L-  . 
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(A26) 


The  solution  to  Eqn.  (A24}  does  not  represent  the  usual  diffusion 
controlled  situation  because  In  order  to  meet  the  symmetry  condition  at  the 
centerline  we  must  give  up  the  physically  sensible  boundary  conditions  at  the 
wall.  Perhaps  variable  coefficients  and  the  Inclusion  of  recombination  can 
remedy  this  problem  with  the  boundary  conditions.  Or  perhaps  we  should  be 
willing  to  give  up  axial  symmetry  and  consider  some  sort  of  contraction 
asymmetry  which  would  likely  have  to  Include  some  reconbinatlon  because  of 
the  higher  pressures.  This,  however  shall  not  be  pursued  In  this  report. 
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APPENDIX  B - INFLUENCE  OF  TWO-BODY  REC(»!BINATION  OF  SCHOTTKY  SOLUTION 

The  diffusion  equation  given  in  the  main  text  (Eqn.  14)  Is  not  suitable 
for  elliptic  funclton  representation,  therefore,  a different  form  Is  developed 
below. 

The  diffusion  equation  with  constant  coefficients  Is  (without  considera- 
tion of  whether  shall  be  used) 

.2  _ 

- D — y ■ V n - on  (Bl) 

* dy^  ^ 

2 

Defining  a stationary  state  v^^n*  - on*  ■ 0 , or 

n*  ■ v^/o  (b2) 

We  can  nondlmenslonallze  the  equation  with 

n ■ n/n*  and  y * y/L 


Obtaining 


d n “ *2 

— y - n - n 

dy^ 


n'  (o)  “ 0 and  n*6,  0<6<1 


n(a/L  ) • 0 
n 


It  Is  clear  that  the  boundary  condition  at  the  centerline  make  up  for  the 
disappearance  of  c In  Eqn.  (B3) . Now  Eqn.  (B3)  can  be  Integrated  once  to 


yield 


I 


FCiPik)  Is  tabulated.  Note  that  for  6 > 1 , we  get  a double  root  (a  - b) 
and  the  equation  Is  not  an  elliptic  Integral  (though  Its  solution  can  still  be 
found  from  the  Tables).  When  S ~ 1 we  get  the  boundary-layer  behavior 
described  by  Eqn.  (19)  or  (20)  In  the  main  text. 

Except  for  6 > 1 , Eqn.  B3  represents  an  eigenvalue  problem.  In  order 
to  check  this,  we  solve  for  the  value  of  the  Integral  at  the  wall  where 
fl  ■ 0 . Thus 


For  very  small  values  of  6 , we  obtain  the  cosine  solution,  namely,  | 

2 1/2  ! 
\ * (Tr/2)  where  L ■ i 

n 

As  E Increases,  we  obtain  progressively  larger  values  until  at  3 > 1 we 

i 

get  X ■ <»  ) . Figure  B1  shows  this  behavior  for  the  eigenvalue.  The  results  | 

shown  In  Fig.  B1  proved  Identical  to  results  from  a computer  program  at  the 
Technlon  (for  the  combined  diffusion  and  heat  dissipation). 

Increasingly  higher  degrees  of  recombination  flatten  out  the  density 
profile  at  the  center  and  Increase  Its  slope  at  the  walls.  The  radial  elec-  | 

trie  field  Is  thus  enhanced.  The  variable  coefficient  solution  Is  not  antici- 
pated to  change  the  above  results  In  any  major  way. 
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